Universality of soft and collinear factors in hard-scattering factorization 
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Universality in QCD factorization of parton densities, fragmentation functions, and soft factors 
is endangered by the process dependence of the directions of Wilson lines in their definitions. We 
find a choice of directions that is consistent with factorization and that gives universality between 
e + e~-annihilation, semi-inclusive deep-inelastic scattering, and the Drell-Yan process. Universality 
is only modified by a time-reversal transformation of the soft function and parton densities between 
Drell-Yan and the other processes, whose only effect is the known reversal of sign for T-odd parton 
densities like the Sivers function. The modifications of the definitions needed to remove rapidity 
divergences with light-like Wilson lines do not affect the results. 
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Introduction Much of the predictive power of QCD 
is provided by universality of the non-perturbative func- 
tions in factorization theorems for hard processes. These 
functions are parton densities, fragmentation functions, 
etc. Whereas the perturbative parts of factorization for- 
mulae can be usefully estimated from first principles by 
weak coupling methods, the non-perturbative functions 
cannot. Universality is the process independence of these 
functions. It allows them to be measured from a limited 
set of reactions, and then used to predict other reactions, 
with the aid of factorization and perturbative calcula- 
tions. 

However, recent developments 0, 0, [jj El show that 
universality is endangered. For example, the Sivers func- 
tion |f| — the transverse single-spin asymmetry of a par- 
ton density — changes sign Q, between the Drell- 
Yan process and deep-inelastic scattering. This is be- 
cause the directions of the Wilson lines necessary for a 
gauge-invariant definition of a parton density depend on 
whether collinear and soft interactions are before or af- 
ter the hard scattering. Important current experimental 
work addresses the associated physics issues. 

Although for parton densities time reversal relates the 
two definitions [l| , the situation is not so clear in general. 
For example, a fragmentation function involves a semi- 
inclusive sum over out-states: 

\A,X,out) (A, X, out|. (1) 

x 

Time-reversal converts these to in-states, and therefore 
does not prove universality with the obvious, process- 
dependent directions for the Wilson lines jj| — although 
the non-universality did not occur in a one-loop model 
calculation 0. Moreover, for hadron production in 
hadron-hadron collisions, Bomhof, Mulders, and Pijlman 
found a jungle of Wilson lines whose universality prop- 
erties are far from clear; see also the comments of Brod- 
sky, Hwang, and Schmidt 0. 



Therefore in this paper we carefully re-examine the 
arguments about Wilson lines to discover the true 
limits of universality, if any. The issues particu- 
larly concern processes that need transverse-momentum- 
dependent (TMD) partonic functions, and we consider 
three such processes: (a) e + e~ annihilation with detected 
almost back-to-back hadrons , (b) semi-inclusive deep- 
inelastic scattering (SIDIS) with measured transverse 
momentum for an outgoing hadron 

USEl. and (c) the 

Drell-Yan process with measured transverse momentum 
[l2j |. For the first process, Collins and Soper showed a 
factorization theorem more than two decades ago . But 
the paper 0] extending the statement of factorization to 
the Drell-Yan process claimed no proof. 

Our new methods show, in addition to the time- 
reversal-modified universality of parton densities, 
that: (1) TMD fragmentation functions are universal be- 
tween e + e~ annihilation and SIDIS. (2) The soft factor 
is universal between all three processes. (3) Universal- 
ity arguments hold even with a redefinition of the non- 
perturbative functions needed to remove the divergences 
due to light-like Wilson lines. Our arguments delimit the 
process-dependent choices of direction that are compat- 
ible with factorization, and for individual processes the 
choice is wider than previously used |J, llll [13J . 

The central technical issue is that a proof of factoriza- 
tion requires an appropriate deformation 14] of momen- 
tum contours out of the "Glauber region" . Allowed di- 
rections of the Wilson lines are those compatible with the 
contour deformation. The possible directions of contour 
deformation are determined by the space-time location of 
soft and collinear interactions relative to the hard scat- 
tering. Hence a careful analysis of one-gluon corrections, 
which forms the bulk of our work in this letter, should 
be sufficient to determine the directions. 

Factorization and Wilson lines Factorization into 
hard, collinear and soft factors results essentially from 
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FIG. 1: Sample virtual one-loop diagrams for e e annihi- 
lation, with fragmentation considered in a spectator model. 



a corresponding set of momentum regions. Individual 
graphs for a cross section have a range of possible lead- 
ing regions of loop-momentum space, and we will follow 



the approach of Collins and Hautmann [13j by using a 
subtractive method that implements a decomposition of 
graphs by possible regions. Full details of the subtractive 
method to all orders have not been worked out explicitly, 
but it should give factorization after a sum over graphs 
and regions. A Ward identity argument is needed to 
disentangle otherwise coupled factors, and it results in 
simple Wilson lines in the gauge-invariant operator defi- 
nitions of the factors. 

Electron-positron annihilation We first consider e + e~ 
annihilation at large Q with two detected hadrons hi 
and in almost back-to-back directions, for which some 
simple graphs with a virtual gluon are shown in Fig. ^ 
(The contour deformation issues arise only for virtual 
gluons.) 
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For the first graph, we make, as in [lflj . a decomposition into four terms: hard, collinear to h%, collinear to h%, and 
soft. Critical for determining the directions of the Wilson lines is the soft term: 

d n l il~l+ 
F+ (P 1 , - (0, r , Or)) F 2 - (P 2 , (1+ , 0, Or) ) ■ 



(27T) 



I 2 + ie 



{-I- + ie) (1+ + ie) {-I- +ie)(l+ - C A l~ + ie) + C B l+ + ie) {1+ + ie) 



MS counterterm. 

(2) 



Here we use light-front coordinates in the center-of-mass with the dominant components of the momenta of h\ and 
h 2 being P x + and P 2 ~ . 



In the hi and h% parts of the graph we have retained 
only the — and + components of the gluon momentum 
I, and we have picked out the dominant components of 
the current to which the gluon couples. We have inserted 
factors of l~ and l + to allow Ward identities to be used 
and then compensated this in the first term in brackets 
by dividing by l~ and l + . With this first term we obtain 
a good approximation to the original graph in the soft 
region, provided that we make a suitable deformation 
out of the Glauber subregion of the soft region. The 
Glauber region is where <C l T . The second and 

third terms are counterterms, to cancel the divergences 
at large positive and negative gluon rapidity that would 
arise if only the first term were used. 

The counterterms have arbitrary positive parameters 
Ca and Cb of order Q 2 /m 2 (with m being a soft scale), 
they are power suppressed in the soft region, and they 
allow a Wilson line definition of the soft factor. The 
arbitrariness of Ca and Cb is exploited by the use of the 
Collins-Soper equation, which controls the C A and 
Cb dependence of the soft and collinear factors, and so 
enables predictions to be made. 

The contour deformation must not cross the final-state 
quark poles in the original graph. We choose to deform 
symmetrically out of the Glauber region: Al + = iw, 



AZ~ = — iw, where w is a suitable positive function of 
the real parts of / M . Compatibility with this deforma- 
tion determines uniquely that the light-like Wilson-line 
denominators are —l~ +ie and l + + ie, and that the signs 
of the Ca and Cb terms relative to their ies are as shown. 
But it does not determine whether the counterterm lines 
are space-like or time-like. A legitimate possibility not 
noticed in is that one or both could be time-like: 
1+ + C A l~ ~ ie, -I' - C B l + - ie. 

There is also an ultra-violet divergence at large It 
which we remove by ordinary renormalization. Neither 
the rapidity divergences nor the UV divergence affect the 
validity of the soft approximation in the soft region. 

In coordinate space, both light-like Wilson lines are 
future pointing, as is intuitively natural. They approxi- 
mate a fast-moving quark and a fast-moving antiquark as 
seen by a slow gluon. (However, the non-light-like lines 
in the counterterms are past-pointing, not so intuitively.) 

We choose space-like lines in the counterterms because 
they are compatible with one of the possibilities for the 
counterterms in SIDIS, and so they allow a proof of uni- 
versality. Furthermore, exact properties of matrix ele- 
ments of Wilson lines are simpler when the gluon fields 
are at space-like separation and hence all commute. 
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Another possible change is to use an asymmetric con- 
tour deformation, such as we choose in SIDIS, i.e., pri- 
marily on /+ only or on l~. But it can be shown that the 
only extra resulting cases for the ie prescriptions violate 
the charge-conjugation relationships between fragmenta- 
tion for quarks and antiquarks. They would therefore 
remove predictive power from factorization. 

In the method of Ji, Ma, and Yuan there are 
no counterterms; instead they use slightly non-light-likc 
lines to cutoff the rapidity divergences. Their square- 
bracket factor would be 



(l+/C A -1- + ie) {-l-/C B + l+ + ie)' w 

The Wilson lines are actually future-pointing, since the 
denominators are on opposite sides of the graph com- 
pared with the corresponding past-pointing space-like 
lines in Eq. J3J). Some differences with our formula- 
tion are inessential power-law corrections. But there are 
also different leading-power contributions to the soft and 
collinear factors. We remark without proof that these 
only occur at large transverse momentum for the gluon 
and therefore amount to a legitimate scheme change. The 
subtractive method provides simpler calculations and a 
cleaner proof of universality. 

Once the soft term is fixed, definite prescriptions for 
the collinear and hard terms follow, just as in |l3l |. so we 
will not present them here. 

There are many other graphs for the process, both with 
different connections of the gluon, as in Fig. ^ and with 
arbitrarily many other lines. In all leading regions, the 
collinear parts are in the final state, so we can continue 
applying the same prescription for the contour deforma- 
tions and for the counterterms. Of course, to use Ward 
identities we must use the same prescription everywhere. 
So we have determined all the Wilson lines. 

Semi-inclusive DIS We now consider similar graphs 
for SIDIS, as in Fig. [21 As concerns the contour defor- 
mation from the Glauber region, the primary difference 
is that in the hi part of the graph the flow of I relative to 
collinear momenta is reversed. So we have denominators 
like (ki+l) 2 — m 2 + ie instead of (fci — I) 2 — m 2 + ie. This 
suggests reversing the contour deformation on l~ , to give 
Al + = iw, Al~ — iw. Then, following [111 ], we might 
reverse the relative signs of l~ and ie compared with Eq. 
(J2J), to get a square-bracket factor 

1 1 

(-/- - ie) (1+ + ie) ~ (-1- - ie) (1+ + C A l~ + ie) 



(-1- -C B l+ -ie) (l++ie)' 

The counterterm lines are now time-like. 

Corresponding time-like lines must also appear in the 
fragmentation function, so that compared with e + e~- 
annihilation, the definitions of both the fragmentation 
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FIG. 2: Sample one-loop diagrams for SIDIS. 




FIG. 3: Sample one-loop diagrams for the Drell-Yan process. 



function and the soft factor are different, and we lose 
manifest universality. Now both of Ca and Cb are large, 
so that we can in fact use space-like lines with our sym- 
metric contour deformation (and similarly in the Ji, Ma 
and Yuan version): 

1 1 

{-I- - ie) {1+ Tie) ~ -ie) (1+ - C A l~ - ie) 



{-I- + C B l+ + te) (l++ie)' 

But even this still differs from the version for e + e~- 
annihilation, so we cannot directly deduce universality. 

Moreover, in the second graph of Fig. [21 the symmetric 
contour deformation is blocked by a trap between initial- 
and final-state poles in the target part of the graph. 
When transverse momenta are of order m, the deforma- 
tion on T is limited to m 2 /Q, not enough to get out of 
the Glauber region. So, as in the proof of factorization 
for diffractive DIS ,15|, we should deform primarily on 
l + , away from the pole(s) in the outgoing struck quark 
and its jet; this typically makes I collinear to the target, 
rather than merely soft. Only small deformations on l~ 
are necessary to avoid the target-side poles both in the 
original graph and in the Wilson-line approximations. 

We can now use exactly the same square-bracket fac- 
tor as in e + e~-annihilation. Since the same direction 
of contour deformation can be applied generally, for the 
Glauber regions for all graphs, the Wilson lines in the 
soft and fragmentation factors are the same as in e + e~- 
annihilation. Thus the soft and fragmentation factors are 
universal between SIDIS and e + e~-annihilation. Space- 
like counterterm denominators, like l + — CaI~ + ie are 
preferred here, since the large positive imaginary part of 
l + assists rather than hinders the deformation of l~ from 
an unphysical pole. 

Drell-Yan process In Fig. 01 we show some graphs for 
the Drell-Yan process. In the first graph, the gluon at- 
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taches to two initial-state lines, so we use a contour de- 
formation opposite to that in e + e~ annihilation. 

But other graphs trap the contour against final-state 
poles in the targ et parts of the graphs. Now to prove 
factorization [16j one can deform away from initial-state 
poles. Crossing target-related final-state poles produces 
extra non-factorizing terms, but these cancel by unitarity, 
after a sum over all hadronic states in the inclusive cross 
section. The argument does not depend on the transverse 
momentum of the lepton pair. 

Therefore we find that we can use initial-state Wilson 
lines for the parton densities and for the soft factor. The 
time reversal argument of applies to all these objects, 
since they all involve only matrix elements of the form 
(tp\AiA2\ip) , where Ai and A 2 are operators and ip labels 
a vacuum or a one-particle state. Such states are the 
same no matter whether they are in- or out-states, so the 
transformation by time-reversal leaves them unaffected. 

This extends the exact universality of parton densities 
and soft factors to the Drell-Yan process, with the excep- 
tion of "T-odd" parton densities, which reverse sign 
as is already known. Our proof now includes Wilson-line 
factors that implement |17| the cancellation of rapidity 
divergences. 

Conclusions and discussion We have shown univer- 
sality of fragmentation functions, soft factors, and par- 
ton densities between e + e~-annihilation, semi-inclusive 
deep-inelastic scattering, and the Drell-Yan process. This 
applies both to the basic definitions with light-like Wil- 
son lines and to the correct definitions with removal of 
rapidity divergences. Regulator lines should be space- 
like. In the Drell-Yan process the lines are reversed com- 
pared with the other processes we considered, but time- 
reversal relates them to the functions for other processes, 
with the usual reversal of sign for the Sivers function and 
other "T-odd" parton densities. 

This eliminates missing elements in the Collins-Soper- 
Sterman formalism ^3] f° r the Drell-Yan process. 

The method of Ji, Ma, and Yuan [llj uses non-light- 
like Wilson lines instead of counterterms for removing 
rapidity divergences. Within this method we find uni- 
versality if the fragmentation function and the soft fac- 
tor in SIDIS have future-pointing space-like Wilson lines, 
contrary to the choice made by these authors. 

Since we need definitions of TMD parton densities that 
differ from the most obvious ones that use light-front 
quantization in light-fro nt g auge, we agree with the con- 
clusion of Brodsky et al. 18] that parton densities are not 
literally probability densities. However, our reasoning is 
different, and builds on the much earlier work of Collins 
and Soper HOH. 

Our work needs extension to hard hadron-production 
processes in hadron-hadron collisions. It is not obvious 



that the extension will succeed. Factorization for these 
processes is at present an unproved conjecture, at least 
for cases where transverse-momentum-dependent parton 
densities and fragmentation functions are needed. 

Even for conventional hadron-hadron-to-hadron fac- 
torization with ordinary, integrated parton densities, 
there is no proof in the literature which correctly treats 
the Glauber region, as far as we know. The factorization 
proofs of Collins, Soper and Sterman |lfi| and Bodwin 
|2fJ are only for the Drell-Yan process. A critical re- 
examination is needed here. 

This work has been supported in part by the U.S. 
Department of Energy under grant number DE-FG02- 
90ER-40577, and by the Sofia Kovalewskaya Programme 
of the Alexander von Humboldt Foundation. 
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